Noname manuscript No. 

(will be inserted by the editor) 



Temperature phase transition model for the 
DNA-CNTs-based nanotweezers 



Anh D. Phan, Lilia M. Woods and N. A. 
Viet 



Received: date / Accepted: date 

Abstract DNA and Carbon nanotubes (CNTs) have unique physical, me- 
chanical and electronic properties that make them revolutionary materials for 
advances in technology. In state-of-the-art applications, these physical prop- 
erties can be exploited to design a type of bio-nanorobot. In this paper, we 
present the behaviors of DNA-based nanotweezers and show the capabilities 
of controlling the robotic device. The theoretical calculations are based on the 
Peyrard-Bishop model for DNA. Furthermore, the influence of the van der 
Waals force between two CNTs on the opening and closing of nanotweezers is 
studied in comparison with the stretching forces of DNA. 
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1 Introduction 

In the past several years, researchers have made much progress in synthesizing 
new materials and developing fabrication techniques necessary for nanoscaled 
device production. This progress has been particularly important for applica- 
tions utilizing physical systems intended for biological and medical purposes. In 
this regard, biophysical devices at the nanoscale open up novel possibilities for 
diagnostic and therapeutic applications. 

DNA and carbon nanotubes (CNTs) are interesting and important systems 
in nanoscience. They have been the subject of many investigations in the 
past two decades 1,2,3,4,5 . DNA is composed of two long polymer strands 
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organized in a double helical structure, where each strand consists of repeating 
units (nucleotides) [6j. CNTs are quasi-one dimensional cylindrically wrapped 
graphene sheets with properties uniquely defined by theregistry dependence 
of the wrapping given by a chirality index (n,m) 7 . Various applications of 
DNA/CNT complexes have been exploited with potential for biosensors [5], 
DNA transporters [ID], and field effect transistors [IT]. The DNA/CNT is 
a composite with complicated structure with temperature dependent motion 
dynamics. Recently, using molecular dynamics simulations researchers have 
proposed molecular tweezers combining DNA and CNTs [T3] - a device with 
further technological and scientific potential. 

A theoretical model of a geometrical soliton of DNA structure was con- 
structed for the first time by Englander [13 (E model). In this model, one 
of the strands of the DNA is represented as a chain of pendula interacting 
with the another fixed similar strand. The E model explains the existence of 
DNA open state due to nonlinear excitations. In addition, the DNA structure 
and dyanmics has been modeled in terms of the Peyrard-Bishop (PB) model 
PP, which has been succssessful in explaining DNA denaturation transitions, 
pre-melting dynamics, and thermal transport. In the PB model, backbone of 
DNA is described as chains of particles with nearest neighboring potentials. 
However, the models ignore the helicoidal structure of the DNA molecule, the 
context of DNA flexibility, and the properties associated with it. 

CNTs are chemically inert and they interact with other materials via long- 
ranged dispersive forces, such as van der Waals (vdW) forces. The vdW inter- 
actions of graphitic nanostructures can be described via pairwise interatomic 
Lennard- Jones (LJ) potentials [T3]. This approach relies on knowledge of the 
coupling Hamaker constants and it predicts the equilibrium separation cor- 
rectly. The LJ potential has been applied to model mutual interaction between 
various CNTs as well as CNT based devices [151115] ■ 

In this work we investigate the dynamics of hybrid DNA/CNT nanotweez- 
ers by employing the PB and vdW-LJ models. This dynamics of stretching 
in terms of its velocity and acceleartion due to environmental temperature 
changes is investigated. The critical temperature where a melting transition of 
the DNA/CNT takes place is presented. Comparisons between the strength of 
the involved forces showing the temperature-dependent motion is dominated 
by the stretching of the H bonds and bases, while the CNT vdW interaction 
is weaker. 

The rest of the paper is organized as follows: In Sec. II, the theoretical 
structure model, behavior and interactions of DNA-based nanotweezer are 
introduced. In Sec. Ill, numerical results are presented. The conclusions are 
given in Sec. IV. 

2 Model and mathematical background 

The proposed nanotweezer architecture is assembled by attaching the reactive 
ends of two single wall CNTs to the DNA strands as shown in Fig. \T\ The 
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rest of the end C bonds are saturated via H atoms. The size of this hybrid 
is quite large, approximately thousands atoms, thus full quantum mechanical 
atomistic treatment is not possible. The PB model is relatively simple pQ, 
which describes the DNA two strands as a coupled pendulum system. 




Fig. 1 (Color online) Schematics of the DNA/CNTs-based nanotweezer. 



2.1 Model of DNA dynamics 



According to the PB model pQ, the DNA double strand is modeled by two 
parallel chains of nucleotides via nearest-neighbor harmonic oscillator inter- 
actions. The potential for the Hydrogen bonds is also included. The relevant 
Hamiltonian is given as follows [TJE] 



N 



n=l 



-m (u 2 n + i>l) + ^fc {u„ - u„-i) 2 + 



1 2 

-k(v n -v n -i) +V(u n -v n ) 



(1) 



where u n and v n are the nucleotide displacements from equilibrium along 
the direction of the hydrogen bonds for each strand, m is the mass of each 
nucleotide (taken to be the same for each unit), while k is harmonic oscillator 
coupling constant of the nearest-neighbor longitudinal interaction along each 
strand in units of eV/ A 2 . The potential for the Hydrogen bonds between the 
two strands is modeled via a Morse potential V(r) = D[e~ ar — l] 2 . Here, D is 
the dissociation energy and a is a parameter. It is important to note that the 
Morse potential represents the hydrogen bonds between complementary bases, 
the repulsive interactions of the phosphate, and the influence of the solvent 
environment. 
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The dynamics of the system described by Eq.((T]) is conveniently described 
using a set of new variables x n = («„ +v n )/ v2 and y n = [u n — v n )/ v2, repre- 
senting the in-phase and out-of-phase motion of the two strands, respectively. 
Using this separation of variables, the Hamiltonian is decoupled. An important 
point is that y n represents the relative displacements between two nucleoid at 
the site n in different strands. It reflects the stretching of DNA. Here we con- 
sider the out-of-phase displacements stretch of the hydrogen bonds given by 



N 



T^myl + \ k (Vn ~ Vn-if + V (2y n ) 



(2) 



In the case of large number of nucleotides N — > oo and H is independent 
on the particular site n. Perfroming statistical averaging in the canonical en- 
semble, the Schrodinger equation of a single mode y using H y is given by [Tj 

nana 



i a 2 

'2/3 2 kdy 2 



V{2y) <p{y) = e<p(y), 



(3) 



where, /3 — 1/kgT, and fcs is the Boltzmann constant. The exact solution for 
eigenenergies is given [H] 



1 , ( f3k\ 2a D ( 1 



a z ( 1 
Wk{ U+ 2 



(4) 



Eq.(j4]) has a discrete energy spectrum when d — (B/a)VkD > 1/2. This 
allows one to obatin a critical temperature T c — 2\fk~D j {aks), which is con- 
sidered as the melting temperature of DNA. The DNA states are continuous 
for T > T c and discrete for T < T c . For the parameters of DNA, when we 
consider T > 200 K, only the value of n — is taken into account. There is 
no excitation state for DNA in our considerations. 

From this, the ground state eigenfunction and eigenenergy in the thermo- 
dynamics limit of a large system is obtained as [TJ[2] 



(y) = V^2' 



d— 1/2 



^r{2d~\) 



,-V2a 



-(d-l/2)V2ay 



(5) 



1 . {Pk\ a D a 2 

In addition, the system described via Eq. ([3]) can be represented as a quasi- 
particle with a tempereture dependent effective mass to* = h 2 (3 2 k. At room 
temperature, the value of the effective mass is approximately 22.87 too, here 
Too is the rest mass of electron. The average stretching of the hydrogen bonds 
can also be calculated via (y) = J Lp 2 (y)ydy [H2]- 
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The stretching force is determined via the expression 

dV((y)) 



d(y) 



(J) 



To investigate thermal properties of DNA, we heated up and cooled down 
temperature of the bio-systems hollowing an expression T — 1.14i + 300 (K) 
[12] . Here T (K) is the environment temperature, t (ps) is time. Basing on the 
average stretching (y) of the coupling constants pointed out above, the velocity 
v = d (y) /dt and acceleration a = d 2 (y) /dt 2 of the opening of the nanotweez- 
ers obtained by taking the hrst and second derivative of the stretching with 
respect to time, respectively, are presented in Fig. [5] 
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Fig. 2 (Color online) (Color online) The time-dependent velocity and acceleration of the 
opening. 
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For k — 2.10~ 3 eV/A 2 , the velocity of the opening increases and reaches 
to the maximum with the value of 10.34 m/s at around t = 20 s. After that, 
the velocity drops significantly to zero. It refers that the temperature corre- 
sponding to the peak is 322.5 K. On the other hand, initially, the value of 
acceleration is positive and rises to the maximum value 0.187 x 10 12 m/s 2 at 
—5.2 s or 294 K before declining gradually to the negative side, crossing the 
time axis at 19.4 s or 322.5 K, touching the bottom —0.164 x 10 12 m/s 2 at 
around 45 s and continuing to approach to 0. It can be easily explained due to 
the fact that below 322.5 K, the stretching velocity climbs significantly, so the 
acceleration is positive. Zero acceleration, of course, is at the relevant bending 
point of the opening velocity. Above 322.5 K, the unzipping velocity declines 
notably, and is nearly unchanged. Therefore, the acceleration has the negative 
values and goes to zero. 

In the same way, for other values of k — 3.10 -3 eV/A 2 and k = 4.10 -3 
eV/A 2 , the zero acceleration takes place at 88 s and 138 s, respectively. It 
means that the melting temperatures corresponds to 401 K for k = 3.10~ 3 
eV/A 2 and 456 K for k — 4.10 -3 eV/A 2 . As a result, there is a possibility to 
obtain the melting temperature by observing the velocity of stretching. 



2.2 CNT van der Waals interaction 



The vdW interaction between the CNT parts of the DNA nanotweezers is de- 
scribed via the Lennard- Jones (LJ) approximation. This approach is widely 
used in calculating disperssive interactions between graphitic nanostructures 
because of its relative simplicity and satisfactory results in determining their 
equilibrium configurations |14j . The LJ potential is essentially a pairwise ap- 
prximation, and for extended systems, one typically perfroms integration over 
the volumes of the interacting objects. For CNTs, the integration is over the 
surfaces of hollow cylinders with radii corresponding to the radii of the nan- 
otubes. The LJ-vdW potential per unit length for two parallel CNTs with radii 
R\ and R 2 is given by [50] 



V vdW = o 2 J J (~ + dSxdS 2 , (8) 

where A and B are the Hamaker constants corresponding to the attractive 
and repulsive contributions, respectively. For graphitic systems, one typically 
takes the values for graphite A = 15.2 eVA 6 and B = 24 x 10 3 eVA 12 [II) . 
a = 4/v / 3a 2 is the mean surface density of Carbon atoms with a — 2.49 A 
being the lattice constant. Also, the distance between the CNT surfaces is p. 
Perfroming the integration over the length of the two CNTs with radii R\ and 
i?2, the LJ-vdW interaction can be written as [2"U] : 

^Ao 2 R 1 R 2 f 27T f 2 " 1 
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Fig. 3 (Color online) Sketch of van der Waals interaction between two CNTs. 
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where the in-plane distance between two surface elements is defined as r 2 = 
(R — Ri costpi + i?2 cost^) 2 + {Ri sin^i — i? 2 sin(/j 2 ) 2 - The definitions of Rx, 
i?2 , <^2 ! an d r are sketched in Fig. [31 

Then, applying the first derivative with respect to R, we obtain the van 
der Waals interaction force per unit length 

F vdw (R) = -^ L - (io) 



3 Numerical results and discussions 

As a prototype, we take that both CNTs are identical with the chiral vector 
(5,0) and (6,0), and lengths L\ = = 5 nm. The total Hamiltonian for the 
system is composed of two term, that account for the stretching and van der 
Waals interaction - H = H y + V v dw- Because of the relatively weak vdW force 
between the tubes, V vc iw is treated as a perturbation compared to H y . The 
parameters of DNA are D = 0.33 eV and a = 18 nm -1 . It is important to 
note that <po(y) and Eq in the previous section is the wave function and energy 
of the ground state of DNA without the presence of CNTs. 

In Fig. 2] and Fig.[5j we show results for the CNT vdW perturbative force 
correction as a function of tempertaure and the stretching force. Fig. [S] indi- 
cates that F s decreases as T increases. The stretching force goes to zero at 
the critical temperature since the properties of DNA change when T reaches 
toT c . 

Obviously, the wave function is temperature-dependent, so the energy and 
energy shift are functions of temperature. The value of £o for three values of k 
at this range of temperature varies from 220 meV to 280 meV. It means that 
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Fig. 4 (Color online) The first-order energy is caused by the van der Waals interactions 
between two CNTs. 

the influence of the van der Waals interaction on the wave function and the 
energy in the ground state is minor. We can calculate separately the interac- 
tions of DNA and CNTs. An additional point is that the larger temperature 
is, the smaller the first-order pertubation of energy is. A simple reason for this 
problem is that when temperature increases, two DNA strands are opened [T] 
and it leads to a rapid growth of distance between two CNTs. 

It is remarkable that we have studied the van der Waals interaction and the 
pertubation energy between two parallel CNTs. This configuration also is used 
in order to calculate all of the van der Waals interactions below. Nevertheless, 
in actual cases, we have two crossed CNTs. The dispersion interaction in real 
biosystems is weaker than that in the parallel state. Therefore, we can utilize 
the wave function (po(y) in the following calculations without addional terms 
due to the perturbation theory. 

It is clearly seen in Fig. [31 at the critical temperature T c , the stretching 
force vanishes because two strands of DNA are broken for T > T c . The opening 
force of DNA is very large at low temperature. The smaller the temperature 
is, the smaller distance between two strands is. This force decreases when 
increasing temperature since the separation distance is larger and larger. 

These results have aggrement about the range of magnitude force with 
experimental data and previous calculations [2I3I2E]- The increase of k causes 
to the growth of stretching force due to the fact that the binding of DNA rises. 

Lets consider the interaction between two CNTs attached in the ends of 
DNA. There are several types of DNA existing in nature such as B-DNA and 
Z-DNA. Since the diameter of DNA is approximately 2.37 nm for B-DNA and 
1.84 nm for Z-DNA. We assume that the initial distance between two centers 
of CNTs is 1.5 nm. It is important to note that the van der Waals force is 
attractive at this range of distance and the sign of this force should be minus. 
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Fig. 5 (Color online) The unzipping force as a function of temperature . 



The magnitude of van der Waals interaction between two CNTs is presented 
in Fig. |H1 




Distance (A n ) 



Fig. 6 (Color online) The van der Waals forces between two parallel CNTs (5,0) and (6,0) 
as a function of the separation distance between two centers of CNTs. 



For k = 2.1CT 3 eV/A 2 , if T < 277 K, the stretching force is much larger 
than the van der Waals force of CNTs (5,0) and (6,0) at the initial state. 
Therefore, it is easy to control the opening and closing of DNA by cooling down 
or heating up. At low temperature, the contribution of the dispersion force in 
the movement of DNA strands is minor. However, it can rise to significant 



10 



Anh D. Phan, Lilia M. Woods and N. A. Viet 



role when T > 277 K. We can do the same way with k = 3.10 3 eV/A 2 and 
k = 4.10- 3 eV/1 2 . 

Figure [7] shows the forces between CNT (5,0) and different CNTs at the 
certain distances. In order to control the opening and closing of nanotweezers, 
the van der Waals force is weaker than the stretching forces. It is difficult to 
operate the movement of nanotweezers if two CNTs have large radii. 




Radius (A ) 

Fig. 7 (Color online) The van der Waals forces between CNT (5,0) and another CNT. 

When we heat up the biosystem, two ends of DNA are separated by the 
stretching force. At the larger temperature, the unzipping force is much larger 
than the van der Waals interactions, the nanotweezers are opened. The ob- 
tained results agree with the previous simulation study [12] ■ Therefore, in our 
nanorobots, the movements of CNTs can be controlled by changing temper- 
ature. In addition, the van der Waals interaction between two cylinders is 
proportional to the length of tubes. If we want to have the smaller van der 
Waals interaction, it is possible to choose the length 1 nm or 2 nm. Another 
point is that long CNTs are bent because of the van der Waals interaction. 
As a consequence, the length of tubes should not be large in designing the 
bio-nanorobots. 



4 Conclusions 

The use of intelligence, sensing and actuation nanodevices in surgery, medical 
treatments and materials science is a reality which has become a hot topic in 
the biomedical industry and research in recent years. Bio-nanorobots provide 
further advance not only in the nanotechnology, but also efficient approaches 
for disease treatment. Our studies showed the behavior and architecture of the 
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bio-nanotweezers. The temperature dependence of the opening displacements 
of tweezers is presented and gives researchers some principles to understand 
the operation of DNA-based molecular machines and devices. In addition, the 
velocity and acceleration of the opening and closing tweezers as a function of 
time are speculated. The theoretical calculations are easy to understand and 
agree qualitatively with the previous works. Further research on these systems 
can considerably extend interdisciplinary implications for the technology. 
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